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Abstract 



A mean-field plus pairing model for atomic nuclei in the Fe region was studied 
using a finite-temperature quantum Monte-Carlo method. We present results 
for thermodynamical quantities such as the internal energy and the specific 
heat. These results give indications of a phase transition related to the pairing 
amongst nucleons, around temperatures of 0.7 MeV. The influence of the 
residual interaction and of the size of the model space on the nuclear level 
densities is discussed too. 
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I. INTRODUCTION 



Quantum Monte-Carlo methods offer an interesting way to study fermionic many-body 
problems. Accurate calculations of ground-state properties have been performed for light 
nuclei, using Variational and Diffusion Quantum Monte-Carlo methods [|TJ. Also for the 
nuclear shell model, Monte-Carlo methods are very useful |2j. They allow to do calcula- 
tions in much larger model spaces than conventional techniques based on diagonalization. 
Furthermore, they are very useful for the study of finite-temperature properties of atomic 
nuclei |3|||]. The same method can also be used as a starting point to calculate nuclear level 
densities [§,§J. We have used a variant of this method to study thermodynamical properties 
of nuclei in the Fe region with a model based on a mean-field plus pairing Hamiltonian. At 
present, we are limited to this too simple model, because our quantum Monte-Carlo method 
is still in a developing phase and because we are limited by our present computer facilities 
(a Dec Alphastation 255/300MHz Workstation and a 200 MHz Pentium PC). Even though, 
at present, our results are not conform to realistic Hamiltonians, they do say a lot about the 
general physics properties. It is a first step to go beyond the mean field approximation. The 
present approach gives results about the possibility of a phase transition related to pairing 
correlations, about the influence of the residual interaction on the level densities and on a 
number of general features of the nuclear many-body structure at finite temperature. 

The combination of a mean-field potential and the pairing Hamiltonian leads to a Hamil- 
tonian H = H m f + Hp. Though this Hamiltonian looks simple, it already leads to a compli- 
cated many-body problem. An often used technique to tackle this problem, is the Bardeen- 
Cooper-Shrieffer (BCS) theory Q. It leads to equations that can be handled easily in a 
numerical way and it has a clear interpretation in terms of quasiparticles. The disadvantage 
is that it gives only an approximate solution, and leads to many-body states in which the 
number of particles is not fixed. For some systems, exact solutions can be found ||. A 
general, accurate solution for this many-body problem is even at present a topic of intensive 
research PJTOf. We have found that our Quantum Monte-Carlo method is a very useful 



method in order to study the ground-state and finite-temperature properties of the nuclear 
pairing Hamiltonian. 

In Section |TJ we introduce the mean-field plus pairing model that we used for our cal- 
culations. In Section |H] aspects of our Quantum Monte-Carlo method are given, with an 
emphasis on the differences with related methods. In Section |V| we present a number of 
results. The role of the pairing-interaction strength, the size of the model space and the 
number of particles is discussed. An estimate is given for the level densities. The lines that 
connect the data points on the figures in this paper are ment to guide the eye. They do 
not correspond to analytical results nor fitted curves, except for the pure mean-field results 
(G=0). Error limits represent 95% confidence intervals. If no error limits are shown, it 
means that they are smaller than the markers of the data points, unless it is stated that no 
error limits were determined. 
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II. A MEAN-FIELD PLUS PAIRING MODEL FOR THE FE NUCLEI. 



For the mean-field potential, a Woods-Saxon potential U(r) is used |TT|: 

where 

V c = Ze 2 /r, r > R c , 

= ~Ze 2 /(2R c )](3-r 2 /R 2 c ), r < R c , 

R c = r c A x l\ 
f(x) = (1 + e*)- 1 with x = (r - r A 1/3 ) /a, 



h 



2.000 fm 2 . 



(1) 

(2) 

(3) 
(4) 



Here, A is the number of nucleons, Z the number of protons. The other parameters are 



taken as in 12 



tq = 1.25 fm, 
r so = 1.25 fm, 



0.65 fm, V = 53.3 + 27(A - 2Z)/A - 0AZ/A 1/3 MeV, 
= 0.47/m, V ao = 7.5 MeV. 



To calculate the mean field and its eigenf unctions, we use the parameters for the nucleus 
2gFe3o- This mean field is used for all nuclei in this particular mass region. For every set 
of quantum numbers I and j, the Woods-Saxon potential is diagonalized in a basis of the 
lowest 60 harmonic-oscillator eigenfunctions with the appropriate symmetry. In this way, the 
single-particle eigenstates and their energies listed in table |, are obtained. Also a number 
of unbound states (with energy E > 0) are obtained. In fact, the Woods-Saxon potential 
exhibits a continuum of unbound eigenstates. Due to the expansion in a finite number of 
basis functions, discrete unbound energy levels are obtained. These can be seen as a discrete 
approximation to the continuum of unbound states. The lsi, lpi, lm, Ids, Ids and 2s i 

2 2 2 2 2 2 

orbitals are considered to be compeletely filled. They form an inert core for the many-body 
problem. The 1/r, 2p3,2pi and If 5 orbitals constitute the valence shell. 

A simple Hamiltonian that accounts for the short-range correlations induced by the 
residual interaction, is the nuclear pairing Hamiltonian Hp 0, that takes the form 



E G 

t=p,n 



fc,fc'>0 



d\' t d]j t dktdkt, 



(5) 



where the operators d\ t create a particle in the corresponding single-particle eigenstates of 
the mean-field Hamiltonian in the valence shell. The index t indicates proton or neutron 
states, k is the time-reversed state of the state k. The notation k,k! > denotes that the 
summation for k and k' should run over states with angular momentum projection j z > 
only. The interaction strength G t depends on the model space and the system under study. 
For the strength of the pairing interaction we take G = 20MeV/56, following the suggestion 
of Bes and Sorensen |13] to take G = 20Mev/A. The same strength is used for protons and 



neutrons, and for all nuclei in the Fe mass region. 
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III. THE QUANTUM MONTE-CARLO METHOD 



The method we used to study the model is based on the Shell-Model Quantum Monte- 
Carlo method. The basic idea of this method is well explained in reference 0. It amounts 
to a decomposition of the Boltzmann-operator e~@ H in a sum of exponentials of one-body 
operators 



e 



-PH = Y^e~ h ^. (6) 



Because of their one-body nature, the terms e - ' 1 ' 3 ' 7 can be handled easily using small matri- 
ces, with a dimension equal to the number of single-particle states considered in the model 
space. The canonical or grand canonical trace of e _/l/3 CT can be calculated with simple al- 
gebraic operations on these matrices. The sum over the auxiliary fields a is then evaluated 
using Monte-Carlo techniques (in casu the Metropolis algorithm WW). 



In our approach, a number of technical aspects are implemented differently from reference 
The main differences are: 

• In order to arrive at a decomposition of the form || we use the Suzuki- Trotter formula 
[JTJJ to separate the one- and two-body part of the Hamiltonian in the exponents. This 
reduces the leading systematic error term to order /? 3 /iV t 2 , with N t the number of 
inverse-temperature slices. Following the prescriptions from reference M, the leading 
systematic error term would be of the order (3 2 /N t . 

• In order to decompose the exponential of the two-body part of the Hamiltonian, we 
use the discrete Hubbard-Stratonovich, described in [0. Compared to the Hubbard- 
Stratonovich transform described in P, p^7|JT8| , it has the advantage that it leads to 



faster matrix operations and to smaller systematic errors. Because of these fast matrix 
operations, we could use a large number of inverse-temperature slices in order to reduce 
the systematic errors (typically N t = 20/3 for the /p-model space and N t = 40/3 for 
the extended model spaces). 

In order to evaluate the canonical trace of e~ h,3 <T , we use the fast an efficient algorithm 



described in reference [19| instead of the number projection technique described in 

The trial steps for the Metropolis algorithm for the sampling of the discrete aux- 
iliary fields a are generated in the following way: a series of M consecutive inverse- 
temperature intervals are chosen randomly to be updated. Because of the permutation 
properties of the canonical trace, we can always shift these inverse-temperature inter- 
vals to the end of the decomposition. The matrix representation U c of the operator 
e -pK j g ca l cu l a ted up to the (Nt — M) th inverse-temperature slice and stored in com- 
puter memory. A number m between 1 and a maximum number is drawn. Then in 
the last M slices m auxiliary fields are drawn. These are then changed randomly. 
The matrix U a > for the altered configuration a' is constructed out of the stored part 
of U a . The canonical trace of U a > is evaluated. Then the configuration a is accepted 
or rejected according to the Metropolis algorithm ||14|| . This procedure is repeated a 
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number of times (typically 7 times), before a new series of M slices is selected. So 
a complete Markov step consists of 7 local updates. Typical values are M = 80 (or 
M = N t for N t < 80) and 1 < m < 160. This scheme allows to update a large num- 
ber of auxiliary fields simultaneously, while requiring only M matrix multiplications 
per udpate (counting the contribution of one inverse-temperature slice as one matrix). 
Only after a complete Markov step all the Nt matrices in the decomposition have to 
be multiplied. 

Observables are evaluated after every 5 complete Markov steps. The values are not yet 
fully decorrelated at this rate (e.g. autocorrelations between 30% and 60% between 
consecutive values for the energy), but leaving a larger interval would not improve the 
performance, because already at this rate the most time-consuming part is the con- 
struction of the trial configurations for the Markov chain. The value of an observable 
A at inverse temperature (3 = l/(kT) is calculated as follows: 



(4 



( 



PH+eA 



e=0 



with 

E W (A) 

A a 



Y,g A a W a 

^- log hr N (e^- E ) 
Tr,v (e*'") . 



e=0 



(7) 

(8) 
(9) 

(10) 

(11) 
(12) 



Here, Tr^r denotes the canonical trace, i.e. the trace over all iV-particle states. The 
mathematical properties of the trace operator are crucial in going from expression ^ to 
expression |8|. The one-body operator resulting from the decomposition ^ with —f3H 
replaced by —(3H + eA, is denoted by hp^e- For the energy (A = H), this can easily 
be implemented as hp^ e = hp- etC . The derivatives are evaluated as 



fjep) -/(-e ) , , 

2eo ' { } 



with eo a small but finite number (typically eo = 1/2048). This way of evaluating 
observables requires a lot of matrix manipulations because the complete matrix U a 
has to be recalculated from scratch for hp >ITt€0 and hp >a _ eo . However, because most of 
the computing time goes to the construction of the trial moves, this has little impact 
on the overall performance. Furthermore, this way of evaluating the observables leads 
to small statistical errors, because it amounts to an insertion of the operator A at each 
inverse-temperature interval, whereas the procedure described in is based on the 
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insertion of A only at the first interval. A special remark concerns the evaluation of 
the specific heat C. This quantity cannot be calculated as the expectation value of an 
observable. We evaluate C after the Monte-Carlo run as 



C = (3 2 (E W (E 2 ) - E W (E) 
with 

d 



Tin (e h ^) 



(14) 



The observable E 2 corresponds to the square of the Hamiltonian: 



E W (E 2 



Tt n H 2 e 



-I3H 



/Tr w (e-^) . 



(15) 



Just like any other Quantum Monte-Carlo method for fermions |21]]22|1 , this method 
suffers from a sign problem at low temperatures. The value of the weights w a can be- 
come negative. This poses a problem for the Metropolis algorithm, because it requires 
that w a can be interpreted as a probability density. The Metropolis algorithm can still 
be used by applying it to the absolute value \w a \. Then we have to treat the sign s a 
of Wa as an observable. Expression || becomes: 



J2a Aa-SfjlWpl 
£<r' S a \w a >\ 

E\ w \(As) 



E 



\w\ 



(16) 



The problem is that the statistical error on this expression scales as 1/Eui(s) (see 
below). Now, for an even number of protons and an even number of neutrons, one can 
exploit a symmetry between states with j z > and states with j z < 0. This symmetry 
guarantees that E|^|(s) is close to 1 even at low temperature jlBl . In reference |Q it 
was claimed that the canonical trace should be exactly equal to one for even particle 
numbers with interactions that lead to matrix representations of the form 



P Q 

-Q* P* 



(17) 



as is the case with a pure pairing force or a pairing-plus-quadrupole force. This is not 
exactly true however, as can be seen from the following counterexample: let U a be a 
4x4 matrix, P and Q are 2x2 matrices, with Q = and 



P 




(18) 



where e is a real number. Then the grand canonical trace is given by 
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Tr M (U a ) 



1 - 2e 2 e 2 ^ + e 4 e 4/3/ \ (19) 



which is positively definite for any value of the chemical potential \x. The canonical 
2-particle trace is given by the coefficient of e 2 ^ M : 



Tr 2 (U a 



2e 2 , (20) 



which is negative. The average sign for some of the systems we studied is shown in 
figure |l| as a function of the inverse temperature (3. The sign for the even particle 
numbers shows a dip around (3 = 2MeV~ l . This is caused by the fact that for a 
fraction of the U a the canonical trace becomes negative, like in the above example. 
However, the symmetry between j z > states and j z < states guarantees that 
the contribution to the trace of the 'fully accompagnied' states, in the sense defined 
in reference ||, will always be positive. These states dominate the low temperature 
trace for even-even nuclei. Therefore the sign tends to one again at low temperatures. 
For odd particle numbers, E| M |(s) tends to zero at low temperatures. So there the 
statistical errors explode. We find that even for odd nuclei we can do calculations at 
j3 = AMeV -1 , which is enough to cool the system almost completely to its ground 
state. For the model studied here, the average sign posed no problems. However, for 
more complicated Hamiltonians, one probably will have to fall back on extrapolation 
techniques as the ones described in references |2J and [P3| , in order to overcome the 
sign problem. 

Because the Metropolis algorithm leads to correlations amongst successive values for 
the observables, care has to be taken to establish accurate error limits. In order to 
get rid of the correlations, 50 independent Markov chains are run for each calculation. 
These Markov chains typically consist of some 600 thermalization steps and 3000 
sampling steps, with an evaluation of observables every 5 th step. 

So the chains are rather short, but long enough to make sure that the computing time 
is not dominated by the thermalization steps. This leads to 50 independent estimates 
Ax, . . . , A 5Q for the quantity {A) p. To obtain the final estimate, we take the weighted 
average of these values, with the average signs si, . . . , s 50 as weights. So we obtain an 
estimate A for (A)p: 

A= E 1 As 1 ^ (21) 

If the Ai and Sj are obtained with enough precision, then a good estimate for the 
statistical error on [H] can be obtained from the following expression: [23 



var(x/y) = jvar(x) — 2E(x/y)cov(x, y) + E 2 (x/y)var(y)| /E 2 (y), (22) 

(var and cov denote the variance and covariance). Taking x = As and y = s, leads to 
an estimate for the variance on A: 



6 



var(yl) 



(23) 



Under the assumption that A is almost normally distributed, which is a good approx- 
imation because A is a weighted average of 50 independent values, this expression 

allows us to determine a 95%-confidence interval A — 2\/var(A), A + 2\/vai(A) for 



A - 2^/var(A), A + 2^var(A) 

{A) p. Expression |2^ also demonstrates that the statistical error is inversely propor- 
tional to the average sign E\ w \(s). 



IV. RESULTS 
A. Proton and neutron contributions 

Some thermodynamical properties of the pairing model for 2fiF e 30 were studied using the 



Quantum Monte-Carlo method presented in section |HI| . Because the proton and neutron 
systems are not coupled to one another, separate results for both particle types are obtained. 
The internal energy of the total system and the contributions of the proton and neutron 
subsystems are shown as a function of temperature in figure 0. The same is done for 
the specific heat in figure |3|. The neutrons contribute more to the internal energy than 
the protons do, because there are more valence neutrons than valence protons. This also 
leads to a slightly stronger peak in the specific-heat curve for neutrons than for protons. 
Qualitatively, there is no big difference between the thermodynamical properties of both 
subsystems. This is not the case at lower values of the interaction strength G. 



B. Dependence on the pairing interaction strength G 

We have studied the pairing model for 26^30 for several values of the pairing interaction 
strength. Calculations were performed for 10 neutrons in a shell with 20 valence states (l/i, 
2p3,2pi and l/s orbitals) and for 6 protons in the same shell. 

The neutron energy as a function of temperature is shown in figure |j. The energy scale 
is chosen such that the inert core has zero energy. The fact that the energy does not go 
to much higher values as the temperature increases, is due to the limited size of the model 
space: not enough high-lying states are included. As we shall discuss later on, the results for 
T > 1.5MeV are not physical anymore. For larger values of G, the system is more strongly 
bound. Furthermore, when raising the temperature, the system stays in its ground state 
longer than for smaller values of G. This indicates that there is an energy gap between 
the ground state and the first excited state proportional to G, as is expected from BCS 
theory [|7|]. The neutron specific heat, as a function of temperature, is shown in figure |5|. 
With increasing strength G, the peak in the specific heat curve shifts to a slightly higher 
temperature and becomes more pronounced. In general, peaks in the specific heat can be 
interpreted as signs of a phase transition. We see here that the pairing correlations, for 
G > 20MeU/56, seem to induce a phase transition in the system. 
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Analogous calculations were done for protons. The proton energy as a function of tem- 
perature is shown in figure |6|. The proton specific heat is shown as a function of temperature 
in figure |7|. The same discussion as for the neutron results, applies here. There is, however, 
a striking difference in the specific-heat curve for low values of G: a second peak develops 
around T = 0.2MeV for G = 10MeV/56. At this value of the pairing strength, the broad 
peak in the specific-heat curve, around T = 0.8MeV, coincides with the peak in the the 
specific-heat curve for a pure mean field. This peak is related to the condensation of the va- 
lence particles in the lowest energy levels of the valence shell (the If 7 orbital). The smaller 
is entirely due to pair correlations, that develop among the 6 particles in the If 7 orbital. In 

figure H the expectation value of the pairing interaction operator Hp is shown as a function 
of the temperature T. While the system with G = 20MeV /56 reaches full pairing strength 
at temperatures T < OAMeV, the system with G = lOMeV/56 comes to this regime only 
at values of T < 0.2MeV. In figure |9|, the number particles in the If 7 orbital and the 
number of particles in the other orbitals are shown. For the system with G = lOMe V/5Q, it 
is observed that approximately all 6 particles occupy states in the If 7 orbital for values of 
T < 0A5MeV. The fact that the pairing correlations reach their maximum for this system 
only at values of T < 0.2MeV, means that the system passes through two phases as it is 
cooled: first, the 6 valence protons condense into the 1/r orbital. At T ~ 0A5MeV, this 
stage is completed. If the temperature is lowered further, pair correlations among these 
particles can develop. At values of T < 0.2MeV the system is almost completely cooled to 
its ground state. For the system with G = 20MeV/5Q, the occupation of the If 7 orbital 
reaches a maximum of about 5.3. The particles always remain spread over all the valence 
orbitals, because now the pairing interaction is strong enough to scatter them out of the If 7 
orbital, even in the ground state. 



C. Dependence on the size of the model space 

For the description of the high-temperature properties of the system, the model space 
given by the fp shell is too small. At temperatures of a few MeV, valence particles can be 
excited to higher-lying single-particle states, or core particles can be excited into the valence 
orbitals or higher energy states. In order to know up to what temperatures the results that 
we obtained in the fp shell are valid, we performed a number of calculations in larger model 

spaces. First, the 3s 1, 2cta, 2cta and log orbitals are added to the single-particle space. 

222 2 

This leads to a many-body problem of 6 and 10 particles in 42 single-particle states. In a 
second extended model the core states are considered as valence states too. Therefore the 
lsi, 1»3, lm, Ids, and lota, 2s 1 are added. Furthermore, also the 3m, 3m, 4s 1 and I07 

2' ^2' ^2' 2' 2' 2 ' ^$ 2 2 ^2 

orbitals are taken into account. This leads to a many-body problem of 26 and 30 particles 
in 78 single-particle states. 

Because multiple shells are used, the model space of the extended systems contains 
spurious excitations related to center-of-mass motion. Therefore, care has to be taken when 
relating high temperature results to internal excitations of the system. For the largest model 
space (without core), these center-of-mass motions can be interpreted as thermal excitations 
of the collective degrees of freedom. This picture would be physically meaningful in the 
absence of a mean-field potential. The fact that the mean-field potential is localized in 
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space, breaks the translational invariance of the model. Therefore, one cannot separate 
the center-of-mass motion from the intrinsic excitations in a clean way [25|. A consistent 
treatment of spurious states is a topic for further research. 

The results for the energy and the specific heat obtained using these model spaces are 
shown in figure |K] to [13|. If the value for the pairing interaction strength G is not changed, 
then a system with a larger model space will have a lower ground-state energy because the 
larger model space allows stronger pair correlations. In order to obtain a comparable pairing 
energy, a reduced pairing interaction strength of G = lQMeV/56 is used for the extended 
model spaces. For the no-core system, the energy is shifted such that the ground-state 
energy coincides with the ground-state energy of the fp shell system. 

In the largest model space, at high temperatures (T > 2MeV), the specific-heat curve 
coincides with the specific-heat curve for G = 0. In this temperature region, the proton and 
neutron energy are some 5MeV lower than in the G = case. Apart from this shift, the 
energy curves are similar to the G = case. This indicates that, at high temperatures, the 
pairing Hamiltonian enhances the binding energy but has no effect on the internal structure. 

At lower temperatures, the specific heat curve deviates from the curve for G — 0, because 
pairing correlations develop. By comparing the results for the fp shell and the first extended 
model space, we see that even the fp shell is too small to describe the system at temperatures 
T > 1.3MeV. In order to compare with the results for the second extended model space 
around temperatures of lMeV, the pairing interaction strength G ought to be reduced 
somewhat more for the latter model space. The vanishing of pair correlations with increasing 
temperature, starting from T ~ lMeV, was also observed in shell-model quantum Monte- 
Carlo calculations for |gFe28 based on more realistic interactions The interesting topic 



of proton- neutron pairing in N = Z nuclei pq| , could of course not be adressed in our too 
schematic model used at present. 

The schematic form of the interaction also urges for a caveat: in the larger model spaces, 
the pairing interaction scatters particles into Woods-Saxon orbitals that extend into regions 
where the mean field vanishes (the orbitals beyond the core and the fp shell). Because the 
pairing interaction introduces correlations in the angular coordinates but not in the radial 
coordinates, one might expect the ground state wavefunction to extend too far into space. 
A more realistic interaction would localize the wavefunction more closely to the core. Our 
calculations show that these extended orbitals contribute only about 1 % of the ground state 
density. Therefore we think that this has only a minor effect on the energy-like observables 
we are interested in here. 



D. Dependence on the number of particles 

We studied systems with various numbers of neutrons in the fp shell: |gFe28, 26^29, 
2gFe 30 and 26^31 are modelled by considering 8, 9, 10 and 11 neutrons in the fp valence 
shell, respectively. For the systems with 9 and 11 neutrons, the sign problem limits accurate 
calculations to values of T > 0.25MeV. Fortunately, this temperature is low enough to get 
a good approximation of the ground state. The neutron energy E n for the various systems 
is shown as a function of temperature in figure [T3|. The proton internal energy is not shown 
because it is equal for all four systems and it is already given in figure |6|. While at high 
temperature the energy curves are equidistantly spaced, with an interval of about 9MeV, 
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there is a relative shift to lower energies for the systems with 8 and 10 neutrons at low 
energy. This is because the pairing correlations are stronger for the even systems than for 
the odd systems at temperatures below lMeV. The shift in the energy for the system with 
10 neutrons can be quantified as 

AE W = E ^ + En11 - E nW , (24) 
with E n g, E n io, E n n the neutron energies for the systems with 9, 10 and 11 valence neutrons 



respectively. The quantity AEio is shown as a function of temperature in figure |15| The 
ground-state energy shift was calculated analogously to expression |4], with the energies 
replaced by the mass excesses given in reference [B^]. A value of 1.77QMeV is obtained. Our 
quantum Monte-Carlo results approach this value remarkably well at temperatures below 
0.5MeV. 



E. Level densities 

Because the internal energy is related to the derivative of the logarithm of the partition 
function Zp and because Z@ is the Laplace transform of the level density g(E) of excited 
states, the results presented above also give information about the level density. The par- 
tition function can be obtained by numerical integration of the internal energy. Then one 
should apply an inverse laplace transform on Zp. Because of the statistical errors on the in- 
ternal energies and hence on Zg, this is however an ill-posed problem. A good approximation 
(at high enough energies) is given by the saddle-point approximation: 

9(E) = . (25) 

Our quantum Monte-Carlo method gives accurate results for E and j3~ 2 C. So the level den- 



sity in the saddle-point approximation is easily obtained. Figure [16] shows the level densities 



derived from the results shown in figures |10| [II], [12] and [13]. The results for the smallest 
model space only go up to excitation energies of 25 MeV. They are in good agreement with 
the results for the first extended model space, up to energies of 20 MeV. However, the results 
for the largest model space deviate from these, even at energies below 20 MeV. This indicates 
that it is important to consider also core excitations when calculating level densities. Figure 



16| also shows the level-density curve from a backshifted Bethe formula cited in reference 
H, with parameters a = h.&QMeV~ l and A = l.38MeV. This parametrization was fitted 
to experimental data p8j in order to reflect finite-temperature properties at temperatures 
between 10 7 an 10 10 i^. Clearly, the level densities are shited too much to lower energies by 
the residual interaction we considered. Shell-model Monte-Carlo calculations with a more 
realistic interaction, but in a limited model space (the Ifz2p32pi lfslgs shell), result in a 
better agreement with the backshifted Bethe formula at energies between 5 and 20 MeV H . 



Figure 17 compares the level density for the largest model space with the level density 



for the mean-field Hamiltonian. For the mean-field case, the level density can be calculated 



exactly using a Monte-Carlo method devised by Cerf [29]. We performed such a calculation 



for the largest model space in order to compare it with the saddle-point approximation. 
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For this largest model space, the saddle-point results follow closely the exact level density 
at high energies. At lower energies, the exact level density shows a structure with a lot of 
peaks, because of the discrete structure of the spectrum. In the saddle-point approximations 
these peaks are absent. If these peaks are smeared out with a width of 0.5 MeV, the curve 
coincides with the saddle-point approximation even at excitation energies as low as 1 MeV. 
Thus the saddle-point approximation gives a good 'smoothed' estimate. In the smaller model 
spaces the agreement is less good. Figure |T7] shows that the residual interaction shifts level 
densities to lower energies. This shift is largest at low energies, because there the pairing 
correlations are strongest. 



V. CONCLUSION 

We conclude by stating that our Quantum Monte-Carlo method offers a powerful tool 
for the study of the nuclear pairing model. We have put emphasis on the thermodynami- 
cal properties. Occupation numbers and the pairing gap can be calculated too using this 
method. Main advantages over other methods are that many-body correlations are taken 
into account exactly, particle numbers are constant and finite temperature results can be 
obtained. The major disadvantage of the method is that spectroscopic information can 
only be obtained indirectly. Finally, we remark that our calculations indicate that pairing 
correlations are important only at low temperature ( below lMeV) and at low excitation 
energies, though they do enhance the binding energy. A signal of a phase transition related 
to pairing is found at temperatures around 0.7MeV. The pairing interaction also shifts the 
level density towards lower energies. Furthermore our results show that it is necessary to 
work in very large model spaces, that also include core excitations, if one wants to calculate 
accurate level densities. The work presented here constituted a part of the Ph.D. thesis of 



S. Rombouts [30 . 
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TABLES 



orbital 



single-particle energies (MeV) 



protons 



neutrons 



2 

1»3 

2 

lpi 

2 

Ids 

2 

Ida 

2 

2si 

2_ 



-34.7106 
-25.3351 
-24.0715 
-15.0034 
-12.7911 
-12.3511 



-42.0333 
-32.2120 
-31.1979 
-21.5607 
-19.6359 
-19.1840 



l/l 

2 

2 P 3 

2 

2pi 

2 

1/5 



-4.1205 
-2.0360 
-1.2334 
-1.2159 



-10.4576 
-8.4804 
-7.6512 
-7.7025 



3si 
2ds 
2da 

2 

l9£ 

2 

3j»3 

2 

3pi 

2 

4s 1 

2 

I91 



4.7316 
5.6562 
6.1324 
6.6572 
6.6663 
6.7469 
8.9016 
9.1386 



-0.3861 
0.2225 
0.9907 
0.5631 
2.5931 
2.6915 
4.4706 
3.5488 



TABLE I. Single-particle eigenstates of the Woods-Saxon potential with the parameters as 
described in section H 
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FIGURES 
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11 neutrons 



p [MeV 1 ] 

FIG. 1. Average sign s as a function of inverse temperature f3 for 8, 9, 10 and 11 neutrons in 
the l/i2p3 2pil/5 shell, G = 20MeV/56. 
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protons 
neutrons 
protons+neutrons 



0.2 0.4 0.6 0.8 1 1.2 

T [MeV] 

FIG. 2. Internal energy £ as a function of temperature T for 6 protons and 10 neutrons in the 
l/l2p3 2pil/5 shell, G = 20MeW56. The energy scale is adapted such that the total, proton and 

2 2 2 2 

neutron internal energies all tend to at low temperature. 
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FIG. 4. Neutron energy E n as a function of temperature T for 10 neutrons in the I/7 2p± 2pi If 5 

2222 

shell, for various values of the pairing strength G 



o 



G=0 MeV 
G=1 0/56 MeV 
G=1 5/56 MeV 
G=20/56 MeV 
G=25/56 MeV 




0.8 1 1.2 
T [MeV] 

FIG. 5. Neutron specific heat C n as a function of temperature T for 10 neutrons in the 

l/l2p3 2pi I/5 shell, for various values of the pairing strength G 

2222 
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G=0 MeV 
G=1 0/56 MeV 
G=1 5/56 MeV 
G=20/56 MeV 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

T [MeV] 

FIG. 7. Proton specific heat C p as a function of temperature T for 6 protons in the 
l/z2ps2pi If 5 shell, for various values of the pairing strength G 

9 9 O 9 
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G=1 0/56 MeV 
G=20/56 MeV 







1.4 1.6 1.8 



0.2 0.4 0.6 0.8 1 1.2 
T [MeV] 

FIG. 8. The expectation value of the pairing-interaction operator (#2} as a function of tempera- 
ture T for pairing strength G = lOMeF/56 and G = 20MeV/56 and 6 protons in the \fi2p\i2p\ If 5 
shell. No error limits were determined. 



> 

CD 
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C 

V 




IIIIIIIIH-j-+ + 







1.4 1.6 1.8 



0.2 0.4 0.6 0.8 1 1.2 
T [MeV] 

FIG. 9. The number particles in the If 7 orbital (full line) and in the other orbitals (dotted 

line) as a function of the inverse temperature 8, for 6 protons in the Ifz2pa2pilf5 shell. No error 

2222 

limits were determined. 
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-60 




-110 1 1 1 1 1 1 1 1 1 

0.5 1 1.5 2 2.5 3 3.5 4 

T [MeV] 

FIG. 10. Neutron energy E n as a function of temperature T for 10 neutrons in the 
l/l2p3 2pil/5 shell (a), for 10 neutrons in the first extended model space (b) and for 30 neu- 
trons in the second extended model space (c). The dashed line gives the result for the second 
extended model space without pairing {G = 0). The results for the largest modelspace are shifted 
such that curve (b) and (c) coincide at low temperature 
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0.5 1 1.5 2 2.5 3 3.5 4 

T [MeV] 

FIG. 11. Neutron specific heat C n as a function of temperature T for 10 neutrons in the 
l/l2p3 2pi I/5 shell (a), for 10 neutrons in the first extended model space (b) and for 30 neutrons 
in the second extended model space (c). The dashed line gives the result for the second extended 
model without pairing (G = 0). 
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CD 
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T [MeV] 

FIG. 12. Proton energy E v as a function of temperature T for 6 protons in the Ifr2p32pilf5 

2 2 2 2 

shell (a), for 6 protons in the first extended model space (b) and for 26 protons in the second 
extended model space (c). The dashed line gives the result for the second extended model without 
pairing (G = 0). The results for the largest modelspace are shifted such that curve (b) and (c) 
coincide at low temperature 




--] 



0.5 1 1.5 2 2.5 3 3.5 4 

T [MeV] 

FIG. 13. Proton specific heat C p as a function of temperature T for 6 protons in the 

l/z2p3 2pi l/s shell (a), for 6 protons in the first extended model space (b) and for 26 protons in 

2222 

the second extended model space (c). The dashed line gives the result for the second extended 
model without pairing (G = 0). 
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FIG. 14. Neutron energy E n as a function of temperature T for systems with 8, 9, 10 and 11 
neutrons in the Ifi2ps2pilf5 shell, (G = 20MeV/56). 
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FIG. 15. Neutron energy shift AE'io as a function of temperature T for a systems with 10 
neutrons in the Ifi2ps2pilf5 shell, (G = 20MeW56). The dashed line indicates the experimental 

2 2 2 2 

value of the ground-state energy shift for 26^30 • 
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/./ 6+10 nucleons in 20+20 states, G=20/56 MeV 

/ 6+1 nucleons in 42+42 states, G=1 6/56 MeV 

26+30 nucleons in 78+78 states, G=1 6/56 MeV 
1 u > backshifted Bethe Formula 



10 20 30 40 50 

E [MeV] 

FIG. 16. Level densities for several model spaces. The ground-state energy in each model space 
was shifted to MeV. The backshifted Bethe formula was fitted to experimental data for 56 i* 1 e. 
(for temperatures up to 10 10 K). 




-1040 -1020 -1000 -980 -960 

E [MeV] 

FIG. 17. Level densities for the largest model space, with and without pairing interaction. The 
energy scale reflects the sum of the single-particle energies as listed in table | and the additional 
binding energy stemming from the pairing interaction. 
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